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ABSTRACT 



Aims. We performed a spatially resolved spectral X-ray study of the pulsar wind nebula (PWN) in the supernova remnant GO.9+0.1. 
Furthermore, we modeled its nonthermal emission in the X-ray and very high-energy (VHE, E > 100 GeV) y-ray regime. 
Methods. Using Chandra ACIS-S3 data, we investigated the east- west dependence of the spectral properties of GO.9+0.1 by calcu- 
lating hardness ratios. We analyzed the EPIC-MOS and EPIC-pn data of two on-axis observations of the XMM-Newton telescope 
and extracted spectra of four annulus-shaped regions, centered on the region of brightest emission of the source. A radially symmetric 
leptonic model was applied in order to reproduce the observed X-ray emission of the inner part of the PWN. Using the optimized 
model parameter values obtained from the X-ray analysis, we then compared the modeled inverse Compton (IC) radiation with the 
published H.E.S.S. y-ray data. 

Results. The spectral index within the four annuli increases with growing distance to the pulsar, whereas the surface brightness drops. 
With the adopted model we are able to reproduce the characteristics of the X-ray spectra. The model results for the VHE y radiation, 
however, strongly deviate from the H.E.S.S. data. 

Key words. X-rays: individuals: GO.9+0.1; ISM: supernova remnants; ISM: individual objects: GO. 9+0. 1; Radiation mechanisms: 
non-thermal; Methods: numerical 
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The composite supernova remnant (SNR) GO.9+0.1 is located 
less than one degree from the Galactic center (GC). It was dis- 
covered with the Molonglo radio telescope by Kesteven (1968). 
Helfand & Becker (1987) found the overall radio morphology to 
be dominated by a luminous core with a diameter of « 2 ' sur- 
rounded by a fainter, but still detectable shell (diameter as 8 ')■ 

The first firm detection of GO.9+0.1 in X-rays has been 
achieved with the BeppoSAX satellite (Mereghetti et al. 1998). 
Follow-up observations with Chandra (Gaensler et al. 2001) and 
XMM-Newton (Porquet et al. 2003) unambiguously identified 
the core of the SNR as a pulsar wind nebula (PWN). The mor- 
phology of GO.9+0.1 seen in the Chandra data set reveals an 
axial symmetry where the symmetry axis lies at an angle of 
~ 165°, measured counterclockwise from north (Gaensler et al. 
2001). An unresolved source (CXOU J174722.8-280915) is 
found along the symmetry axis. Due to the high sensitivity of 
XMM-Newton, Porquet et al. (2003) were also able to detect 
faint X-ray emission from the shell of the SNR. Comparing dif- 
ferent regions in the PWN, they observed a positive gradient of 
the spectral index from east to west. This is interpreted as fast 
rotation with resulting Doppler boosting of the leptons that con- 
stitute the PWN, requiring an anomalously large extent of the 
termination shock with a shock radius of R$ > 1 pc as suggested 
by Gaensler et al. (2001). 

Dubner et al. (2008) have compared newer radio observa- 
tions with a revised XMM-Newton analysis. These authors 
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show that the X-ray PWN almost fills the size of the ra- 
dio core, which indicates a moderate magnetic field (see e.g. 
Gaensler & Slane 2006). They support the axial symmetry sug- 
gested by Gaensler et al. (2001), but do not see any spectral vari- 
ations in the radio regime, contrary to the ones in X-rays reported 
by Porquet et al. (2003). The authors did not find a radio coun- 
terpart of CXOU J174722.8-280915. 

GO.9+0.1 has been detected with H.E.S.S. (Aharonian et al. 
2005) in very high-energy (VHE) y-ray s above 200 GeV. With 
a flux of only 2 % compared to the one of the Crab Nebula, 
GO.9+0.1 was the faintest known VHE y-ray source at that time. 
However, due to the large field of view of H.E.S.S., the tele- 
scopes simultaneously observed GO.9+0.1 when pointing to- 
wards the GC, resulting in a high significance of 13cr due to 
the long exposure. For H.E.S.S., the source appears point-like, 
with an upper limit on the intrinsic angular extent between 1.3 ' 
assuming a Gaussian emission region and 2.2 ' for the emission 
from a thin shell (both at 95 % confidence level). 

In 2009, Camilo et al. discovered PSR J 1747-2809, the pul- 
sar powering GO. 9+0. 1, using the NRAO Green Bank Telescope 
at a frequency of 2 GHz. The authors report a strongly scat- 
tered and dispersed signal with a pulsation period of P — 52 ms. 
The characteristic age of the pulsar (Manchester & Taylor 1977) 
was derived as t c — 5.3 kyr and the spin-down luminosity as 
E — 4.3 x 10 37 erg/s. Notably, this corresponds to the second- 
highest known spin-down luminosity of Galactic pulsars, only 
surpassed by the Crab pulsar. Camilo et al. (2009) estimate the 
true age of PSR J 1747-2809 at 2 - 3 kyr, hence well below the 
characteristic age. The high scattering of the signal increases 
the positional uncertainty to about 3 ', neither proving nor dis- 
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Table 1. Analyzed GO.9+0. 1 XMM-Newton observations. 



Observation 


Year 


Instrument 


Exposures 


(ks) 


ID 




EPIC 


performed* 


net (2 > 


0112970201 


2000 
2000 


pn 
MOS1/2 


11.7 
17.2 


11.3 
17.2 




2003 


pn 


43.7 


26.9 


0144220101 


2003 


MOS1 


49.4 


44.1 




2003 


MOS2 


49.4 


45.3 



Before background screening 
1 After background screening 



proving that CXOU J174722. 8-280915 is the X-ray counterpart 
of PSR J1747-2809. Due to the high absorption and hard spec- 
trum, Gaensler et al. (2001) ruled out the possibility that CXOU 
J174722. 8-280915 corresponds to a foreground star. The au- 
thors estimate that the probability of its emission originating in 
a background active galactic nucleus is « 10~ 2 . 

Up to the present, the distance to GO.9+0. 1 is not very well 
constrained. Many authors (for example, Aharonian et al. 2005; 
Dubneretal. 2008) adopt d = 8.5 kpc, which would imply 
GO.9+0. 1 lying in the direct vicinity of the GC. Applying the 
NE2001 electron model (Cordes &Lazio 2002) to the disper- 
sion measure of PSR J1747-2809, Camilo et al. (2009) obtain 
d ~ 13 kpc, but admit that this model could be substantially in 
error for sources toward the inner Galactic regions. Due to this 
uncertainty, they propose an approximate distance of 8 - 16 kpc, 
leaving this problem unsolved. 

In the subsequent section, we describe the analysis of the 
XMM-Newton observations used for this work. The following 
part contains a detailed study of the spectral properties of the 
PWN in GO.9+0. 1 in X-rays, using data from Chandra, as well 
as from XMM-Newton. In Section 4 we outline the radially sym- 
metric leptonic model applied to this source. The results of the 
modeling for two different lepton injection spectra are presented 
and compared in Section 5. Using the parameters optimized to fit 
the X-ray data, we calculate the inverse Compton (IC) emission 
and compare it with the H.E.S.S. results in Section 6. 

2. XMM-Newton observations and data analysis 

GO.9+0. 1 was observed on-axis twice by the XMM-Newton 
telescope (Jansenetal. 2001) in 2000 and 2003. Table 1 lists 
the properties of these two observations used in our analysis. 
There is also a third data set available (Observation ID (OID): 
0205240101), but the high off-axis position of approximately 
12 ' makes a spatially resolved spectral analysis difficult due to 
the strong distortion of the point spread function (PSF) at large 
off-axis angles. The earlier data set (OID: 0112970201) has al- 
ready been extensively analyzed by Porquet et al. (2003). 

During both on-axis observations, the medium filter was 
used for the EPIC-MOS (Turner et al. 2001) and EPIC-pn 
(Striider et al. 2001) cameras. The MOS cameras were operated 
in the standard full-frame mode in both cases, whereas the pn 
camera was operated in the extended full-frame mode (2000) 
and in the large window mode (2003). 

We analyzed the data using version 9.0.0 of the Science 
Analysis System (SAS) provided by the XMM-Newton Science 
Operations Centre, together with tools from the FTOOLS pack- 
age (Blackburn 1995). To reduce the influence of soft proton 
flaring, we selected good time intervals by applying a maximum 




Fig. 1. Chandra count map in the energy range 3-8 keV. Shown 
in red are the two annulus halves chosen for the hardness ratio 
calculation (centered on the region of brightest emission in this 
case). The point source CXOU J174722.8-280915 is encom- 
passed by a dashed green circle. 



count threshold to the high-energy lightcurve (7-15 keV) ob- 
tained with the standard analysis chain. Values of eight, respec- 
tively two, background counts per second were chosen for the 
pn and for the two MOS cameras. The resulting net exposures 
are also shown in Table 1 . For each data set, we only selected 
good (FLAG = 0) single and multiple (PATTERN < 4 for pn 
and < 12 for MOS) events. Since the source is extended and has 
a moderate surface brightness, it is not necessary to correct for 
pile-up. 

For each camera and observation, we defined an infield back- 
ground region located near the source and on the same CCD 
chip. For extended sources with low surface brightness like 
GO.9+0. 1, a proper treatment of all background components is 
particularly important. Therefore we used the method for par- 
ticle background subtraction as described by Majerowicz et al. 
(2002) which was applied to the infield background, as well as 
to the respective source region. To take the vignetting of the 
telescopes into account, we adopted the weighting method in- 
troduced by Arnaud et al. (2001). 

Channels were grouped to a signal-to-noise ratio of five for 
all spectra. 

3. Spectral properties 

Due to its small point spread function and low background level, 
Chandra is especially well suited to investigating spectral differ- 
ences on a small morphological scale. Using the 35 ks Chandra 
data set (OID: 1036) previously analyzed by Gaensler et al. 
(2001), we calculated the hardness ratios using the counts in- 
side two opposite annulus halves with inner and outer radii 5 " 
and 45 ", respectively. The halves are separated by the aforemen- 
tioned symmetry axis, as illustrated in Fig. 1. They were cen- 
tered once on CXOU J174722. 8-280915 and once on the region 
of brightest emission. We calculated the fractional difference us- 
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Table 2. Hardness ratios of the estimated source count inten- 
sities within the annulus halves. The errors are given at 68 % 
confidence level. 



Center 


Orientation 


HR 


CXOU J174722.8-280915 


East 
West 


O 37+0.04 
"""-0.05 

-0 25 +0 04 

"""-0.04 


Brightest emission 


East 
West 


-0 36 +0 04 

"•-' o -0.04 
A 9Q+0.06 
^■^-0.06 




ing the Bayesian Estimation of Hardness Ratios from Park et al. 
(2006): 

Au — As 

HR= — -, (1) 

An + M 

where As and An are the expected soft and hard source count 
intensities, calculated from the measured ones by taking the dif- 
ferent effective and extraction areas into account. The low and 
high energy bands are 3-5 keV and 5-8 keV, respectively. The 
obtained values are given in Table 2. When centering on CXOU 
J174722. 8-280915, the resulting hardness ratios of the eastern 
and western halves are inconsistent. The results are, however, 
comparable within the statistical errors when the halves are cen- 
tered on the region of brightest emission. 

For the spectral modeling of the XMM-Newton data, we 
used version 12.7.0 of the XSPEC tool (Arnaud 1996). The ex- 
tracted spectra of the pn and MOS cameras of the 2000 and 2003 
observations were fitted in parallel with an absorbed power-law 
model. We used the tbabs absorption model, along with the 
abundances of the interstellar medium from Wilms et al. (2000). 

To obtain a statistically significant result for the absorption 
column density, we first fitted the spectrum of a circular region 
with a radius of 45 ", in the energy range 0.2 - 10 keV. The cir- 
cle was centered on the region of brightest emission. The result 
is N H = (2.25 ± 0.15) x 10 23 cm" 2 . Due to the different abun- 
dances used in our fit, this value is considerably higher than the 
one obtained by Porquet et al. (2003, N H = (1.39 + 0.13) x 10 23 
cm -2 ; both values at 90% confidence level). When using the 
same abundances, the values are equal. 

For the modeling of the emission (see Section 4), we ex- 
tracted spectra of four annuli centered on the region of brightest 
emission. They are illustrated in Fig. 2, which shows a smoothed 
XMM-Newton count map of GO.9+0. 1, merged from the MOS 
and pn data of the 2000 and 2003 observations. This map was 
generated using the SAOImage ds9 tool (Joye & Mandel 2003). 
The properties of the annuli are given in Table 3. They are num- 
bered with increasing distance to the center. Since the appear- 
ance of GO.9+0. 1 in X-rays obviously deviates from radial sym- 
metry at least at higher angular distances, we refrained from ex- 
tracting spectra with greater radii. To increase the statistical sig- 
nificance, we fixed Nh to 2.25 x 10 23 cirT 2 for the individual 
annuli. The results obtained from the fitting, as well as the de- 
rived surface brightness, are also given in Table 3. According to 
the reduced x 1 values, a power-law fit yields a very good esti- 
mation for each of the spectra. As an example, Fig. 3 shows the 
spectrum of the second annulus. The strong absorption below 
~ 3 keV is clearly visible. 

4. The model 

In the following, we model the observed emission of the PWN in 
GO.9+0. 1 with a leptonic scenario. To comprehend the softening 




Fig. 2. Smoothed XMM-Newton count map merged from the 
MOS and pn data of the 2000 and 2003 observations (see 
Table 1) in the energy range 3 - lOkeV. The annuli are encom- 
passed by dashed cyan lines (properties listed in Table 3). 




Energy (keV) 

Fig. 3. Energy spectrum extracted from the second annulus. The 
dark blue and black curves correspond to the pn camera of the 
2000 and 2003 observations, respectively. The four lower curves 
represent MOS 1/2 of both observations. The data were fitted in 
parallel with an absorbed power-law model (tbabs*cflux*po). 



of the spectra with greater distance to the pulsar, we adopt a 
spatially resolved instead of a one-zone model. 



4. 1 . Lepton injection spectrum 

A pulsar with spin-down luminosity E loses part of its energy 
in the form of a lepton-dominated particle wind, called the pul- 
sar wind. The leptons of the pulsar wind are accelerated at the 
wind termination shock and injected into the PWN. The gen- 
eral shape of the injected lepton spectrum is commonly as- 
sumed to be a broken power law with the following form (e.g. 
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Table 3. Results obtained by fitting a power-law spectrum to the XMM-Newton data of GO.9+0. 1. All errors are quoted at the 1 cr 
level. 



Annulus 


Radius (") 


r 


Unabsorbed flux'*' 


Surface brightness**' 


^/dof 


No. 


inner 


outer 




(10~ 12 ergcnr 2 s" 1 ) 


(10~ 16 erg cirr 2 s~'arcsec~ 2 ) 




1 


4 


10 


1.37 ±0.09 


0.63 ± 0.02 


24.0 ± 0.6 


79/71 


2 


10 


20 


1.59 ±0.07 


1.07 ±0.02 


11.3 ±0.2 


97/103 


3 


20 


30 


2.01 ± 0.08 


0.91 ±0.02 


6.8 ±0.1 


87/81 


4 


30 


45 


2.18 ±0.08 


1.18 ±0.02 


3.5 ±0.1 


86/83 



*' Energy range 2.5 - 10keV 



Kennel & Coroniti 1984b; Reynolds & Chevalier 1984): 



\ E b,V 



Go,p 



£b,P 



for £ e < E b p, 



for E e > £ b-P , 



(2) 



with the lepton energy E e , the break energy E b p, the spectral 
indices p\ and p2, and the normalization <2o,p of the spectrum. 
Since the model focuses on the inner part of the PWN close to 
the pulsar where mainly young leptons are expected to contribute 
to the nonthermal emission, we assume that the time scale of 
our modeling is smaller than any potential variability time scale 
of the injection spectrum. Furthermore, the applied model only 
aims to reproduce the measured X- and VHE y radiation, which 
allows the view to be restricted to the part of the lepton spectrum 
where E e > £b,p- 

Spitkovsky (2008) suggests an alternative shape for the in- 
jected lepton spectrum using two-dimensional particle-in-cell 
simulations. The best-fit solution to the shape of the lepton pop- 
ulation downstream of the termination shock is given by 



'go.si^e exp(-^-) 



Qs(E e ) = 



for E e < £ bjS , 



e ,si£eexp(-^ 7 ) + e ,s2£; as 

for Eh S < E e < E cut , 



go.si^eexp^-^-) + Q ,s2E e " s exp(- 
for E e > E cut , 







(3) 

with AEi the width of the relativistic Maxwellian that defines the 
low-energy domain of the spectrum, as the index of the power 
law, E cut and A£ cut its cutoff energy and width, and Qo,s\ and 
Qo,S2 the normalizations of the respective functions. According 
to Spitkovsky (2008), as lies in the range 2.3 - 2.5. The ratios of 
the relevant quantities obtained from the fit are E\, t s/AE\ ~ 7, 
£cut/£b,s ~ 7.5, and E cu t/AE cu t ~ 3. For our modeling, we 
used the same ratios to preserve the shape of the spectrum. 
Furthermore, we fixed as to a value of 2.4. 

In the remainder of this paper, we refer to the case of the 
power-law injection spectrum as Scenario I, whereas Scenario II 
corresponds to the injection spectrum from Spitkovsky (2008). 
The other parts of the model are identical for both scenarios and 
are laid out in the following. 

The total amount of spin-down power transferred into the 
energy of the modeled leptons is 



r 



Qp/s(E e )E e d£ e = rjE, 



(4) 



where 77 denotes the conversion efficiency of the pulsar. Using 
this equation, the normalization of the lepton spectrum can 
be calculated when the lower and upper integration limits are 
known. We chose £e,min = 1 erg, which is well above typical 
break energies of the broken power law (£ b ,p); however, this 
value is still low enough not to affect the synchrotron or IC 
emission in the observed energy ranges, as can be derived from 
de Jager & Djannati-Atai (2009). Independent of the accelera- 
tion mechanism, there are two limits that constrain the maximum 
energy to which leptons can be accelerated at the termination 
shock (de Jager & Djannati-Atai' 2009). The first one states that 
the gyroradius of the charged particles constituting the pulsar 
wind must be smaller than the shock radius R$. When defining 
e = Rl/Rs with R^ the gyroradius, the maximum energy of the 
leptons is given by (de Jager & Djannati-Atai' 2009) 



(T 



\+cr c 



(5) 



In this equation, k denotes the magnetic compression ratio at the 
shock and cr the magnetization parameter corresponding to the 
ratio of the magnetic and the particle energy outflow. For want 
of better estimates and due to the partly strong correlation of 
the parameters, we fixed e = 1 . According to the second limit, 
a charged particle reaches its maximum energy when the syn- 
chrotron losses become as strong as the energy gain. This state- 
ment can be rewritten to (de Jager & Djannati-Atai 2009) 



43.7 B~f erg, 



(6) 



with Z?s,g the magnetic field strength at the shock in units of 
G. The lower value of both constraints is then used as £ e up for 
the power-law shape or as E cut when using the spectrum from 
Spitkovsky (2008). In the latter case, £e,up is set to infinity. 



4.2. Outward propagation of the leptons 

The radially symmetric model applied to GO.9+0. 1 resembles 
the one introduced by Schock et al. (2010). The lepton plasma is 
assumed to propagate outwards with a bulk velocity of 



v(r) = v s (y) , 



(7) 



where a is the index of the adopted power law and \'s the ve- 
locity at the wind termination shock. When combining k and 
cr by defining £ = k V°"/(l + °"X the relation for the magnetic 
field strength at the shock is given by (Kennel & Coroniti 1984a; 
Sefako&de Jager 2003): 



(8) 
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Assuming a toroidal magnetic field whose outward propagation 
is directly connected to that of the leptons, the ideal magnetohy- 
drodynamic (MHD) limit on the assumption of a static system 
yields (Kennel & Coroniti 1984a) 



Bvr = BsvsRs = const. 



(9) 



As the leptons propagate outwards, they lose energy, leading 
to a change in the spectral shape. Two fundamental energy 
loss mechanisms have to be considered: synchrotron radiation 
of the leptons and adiabatic energy losses. They are given by 
(de Jager & Harding 1992) 



d£e 
df 



V- Vj.0) - 2.368 x KT 3 (B£ e ) 
3 s 



2 er § 



(10) 



In this equation, the first term corresponds to adiabatic energy 
losses and can be calculated using Eq. 7. The second term de- 
notes the synchrotron losses of the leptons. 

4.3. Implementation and photon emission 

For the numerical implementation of the model, we divided the 
observed part of the PWN into a large number of concentric sub- 
shells in order to simulate the continuous case, exactly as per- 
formed by Schocket al. (2010). Using the previously described 
relations, the lepton injection spectrum is then propagated out- 
wards from one subshell to the next. 

The emitted synchrotron and IC radiation of the corre- 
sponding lepton population can be calculated for each subshell. 
For this, we used the equations given by Blumenthal & Gould 
(1970). Regarding the IC process, essentially three seed-photon 
fields are relevant in the case of GO.9+0. 1: the CMB compo- 
nent, IR photons emitted from local dust, and the starlight com- 
ponent. The CMB spectrum is described well by a blackbody 
distribution for arbitrary locations. As an approximation of the 
IR and starlight components, we used the interstellar radiation 
fields of Porter & et al. (2005) developed for the GALPROP 
code (Strong et al. 2000). The emission of the subshells was 
summed up to shells with the same inner and outer radii as cho- 
sen for the annuli of the X-ray analysis. 




20 25 30 

Distance to Pulsar [arcsec] 

Fig. 4. Evolution of the spectral index of the X-ray emission with 
increasing angular distance to the pulsar. The blue and red points 
correspond to the results obtained with the optimized parameters 
of Scenario I and Scenario II, respectively. The XMM-Newton 
data are shown as black crosses. 



This means that, for example, a superposition of 74 % of the sec- 
ond, 33 % of the third, and 12 % of the fourth shell generates the 
emission of the second annulus. We added up the modeled emis- 
sion from the shells accordingly in order to obtain the resulting 
spectra of the annuli. 



4.5. Parameter optimization 

The previously described model allows us to calculate the syn- 
chrotron emission of the annuli for given parameters p (in the 
power-law case), R$, v$, r\, and a. Afterwards the unab- 
sorbed energy flux is computed in six energy bins of equal width 
(1.25 keV) between 2.5 and lOkeV for each annulus. The re- 
sults are then compared with the measured XMM-Newton data 
by calculating an X 2 value as 



4.4. Projection effect 

Since the model is three-dimensional, a shell actually corre- 
sponds to a hollow sphere in contrast to the annuli of the X-ray 
analysis. Such an annulus is merely the two-dimensional projec- 
tion of the three-dimensional shells. This has two major impli- 
cations. First, only part of the volume of a shell is visible in the 
projection represented by an annulus. Second, outer shells add 
up to an annulus emission. The overall emission from annulus i 
is given by 

v-i dN 



dN 
dE 



/,Ann 



(ID 



with dN/dE\j the spectrum of shell j. This equation includes the 
assumption that the radiation is emitted isotropically and does 
not undergo any absorption inside the PWN. The entries of p,j 
can be written in a matrix of size « x «. For the chosen extraction 
regions from Table 3, n — 4 and the matrix can be calculated to 



V- = 



0.82 0.22 0.07 0.03^ 
0.74 0.33 0.12 
0.59 0.26 
0.59j 



(12) 



X 1 



s 



- F u 



AFn 



(13) 



where F,,xmm and F, moc j denote the measured and modeled flux 
and AFi xmm the statistical error in the corresponding energy bin. 
The parameters of the model are scanned over the allowed range, 
searching for a minimum of X 2 . Splitting the spectra into several 
bins allows optimization of the parameters by taking the overall 
energy flux into account along with the spectral shape. 



5. Results of the modeling 

The parameter optimization was carried out separately for 
the aforementioned injection spectra (see Sect. 4.1). We as- 
sumed a distance of d — 13kpc for both scenarios, following 
Camilo et al. (2009). The upper limit on the radius of the termi- 
nation shock was set to = 4 ". 

Since our model only aims to reproduce the inner part of 
the PWN, we carried out the parameter optimization while ig- 
noring the X 2 values of the fourth annulus. The measured and 
optimized modeled values of the spectral index are presented in 
Fig. 4. Although the modeled index for Scenario I also increases 
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Fig. 5. Surface brightness plotted over the distance to the pulsar. 
The energy range is 2.5 - lOkeV and the blue and red points 
correspond to Scenarios I and //, respectively. The black crosses 
denote the measured XMM-Newton data. 




20 25 30 

Distance to Pulsar [arcsec] 



Fig. 6. Radial dependence of the magnetic field strength for two 
parameter sets of Scenario II with a low X 2 value. The shock 
magnetic field strength Z?s is indicated with a dashed line for 
both cases. 



Table 4. Lower and upper values of the parameters with an X 2 
maximum of 15 % above the lowest obtained value for the opti- 
mization of the first three annuli. 



Parameter 


Scenario I 


Scenario II 


P 


1.7 






0.08 


0.22 - 0.26 


a 


1.00- 1.05 


1.10- 1.15 


>Ps ["] 


2.5-4.0 


4.0 


vs/c 


0.13-0.31 


0.3 


n 


0.4-0.7 


0.08 


Rs [pc] 


0.16-0.25 


0.25 


fis IjuG] 


3.9 - 6.2 


10.7- 12.7 



with growing distance to the pulsar, the slope is too shallow. In 
Scenario II, by contrast, the spectral index exhibits a stronger 
steepening and fits the measured data. The evolution of the sur- 
face brightness is illustrated in Fig. 5. The results of the parame- 
ter optimization for both scenarios agree with the measured data. 
As the parameters of the model are partly correlated, one can find 
several very different parameter sets with similar X 2 values. It is 
thus not possible to designate the fit minimum to a definite pa- 
rameter set. The lower and upper values of the parameters with 
an X 2 maximum of 15 % above the lowest value for the opti- 
mization of the first three annuli are shown in Table 4. However, 
only certain parameter combinations lead to a good fit result. 
The best fit results for the conversion efficiency r\ of Scenario I 
are higher than the ones of Scenario II, but they lie well within 
the physical range below 1 for both scenarios. The velocity at 
the termination shock vs is barely constrained for Scenario I, 
but strongly tends to the theoretically expected value of c/3 
(see e.g. Kennel & Coroniti (1984a)) for Scenario II. (p$ remains 
quite unconstrained for the power-law case, whereas it reaches 
the aforementioned upper limit of 4 " when using the injection 
spectrum of Spitkovsky (2008). The results for this parameter 
should, however, not be overemphasized due to its strong cor- 
relation with vs. When combining Eqs. 7 and 9, the radial de- 
pendence of the magnetic field strength can be calculated for a 
given parameter set. This was performed for two sets of Scenario 



II that both exhibit comparably low X 2 values but different shock 
magnetic fields, as shown in Fig. 6. The values of the shock mag- 
netic field strength and velocity index are B$i = 12.7/iG and 
a\ = 1.15 for the first and Z?s2 = 10.7 /xG and ai = 1.10 for 
the second parameter set, respectively. Since the velocity index 
or is greater than 1 in both cases, the magnetic field increases 
with growing distance to the pulsar. However, we are only fo- 
cusing on the inner part of the PWN. With even greater distance, 
the ideal MHD limit no longer holds, implying that the magnetic 
field is not frozen into the particle propagation any longer and 
might lead to a negative slope farther out. The average magnetic 
field strength is still rather low for both scenarios, in agreement 
with the interpretation of Dubner et al. (2008). Notably the val- 
ues for Scenario II are also in very good agreement with the ones 
found by Tanaka & Takahara (2011, B = 12 - 15/iG). 

6. Implications for the VHE y-ray emission 

As for the synchrotron emission, the IC radiation of the whole 
modeled area can be calculated for a given parameter set. As 
already mentioned in Sect. 1, GO.9+0. 1 was also detected in 
the VHE y-ray regime with H.E.S.S., enabling a comparison 
of the modeled with the measured data. The spectral energy 
distribution with the modeled synchrotron and IC emission for 
Scenario II is shown in Fig. 7. It also comprises the XMM- 
Newton data of the same area, as well as the H.E.S.S. data for 
the whole source (Aharonian et al. 2005). While the synchrotron 
emission matches the measured X-ray data well, the IC contri- 
bution strongly deviates from the H.E.S.S. data concerning the 
spectral shape, as well as the flux, implying that the model is 
not well suited to reproducing the VHE y-ray emission of this 
source. This deviation may come from a number of superimpos- 
ing effects that are discussed in the following. 

Due to the large PSF of H.E.S.S. (* 6', Aharonian et al. 
2006), GO.9+0. 1 appears point-like in VHE y-rays. This implies 
that the modeled and measured VHE y radiation likely originates 
in regions of different size. Assuming a spatial emission distribu- 
tion with p oc exp(-8 2 /2cr 2 om . ce ), Aharonian et al. (2005) obtain 
an upper limit on the source extent of <x soulce < 1 .3 ' at 95 % con- 
fidence level. The IC emission within 30 " (i.e. 0.5 ') around the 
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Fig. 7. Modeled spectral energy distribution of GO.9+0. 1 show- 
ing the synchrotron and IC emission of the inner 30" for 
Scenario II. Also included are the measured XMM-Newton data 
of the modeled area, as well as the H.E.S.S. data of the whole 
source (Aharonian et al. 2005). 

center, corresponding to the modeled part of the PWN, would 
contribute only 7 % to the total VHE y radiation of GO.9+0.1 
in that case. Furthermore, we only modeled the leptons that are 
freshly injected at the wind termination shock of the PWN and 
are confined in the MHD outflow. The lifetime of a lepton scat- 
tering photons to the VHE range, however, is longer than the 
one of a lepton, which emits X-rays in the keV energy range 
(de Jager & Djannati-Atal 2009). Thus a part of the measured 
y-rays is expected to originate in older leptons that are not en- 
ergetic enough to emit X-rays any longer. In addition to that, 
it is also possible that the PWN actually contains a significant 
amount of energetic hadrons in addition to the leptons. If target 
material is present, they would produce 7r°s, which in turn decay 
into two y-rays each, increasing the overall VHE y flux of the 
source. Such a lepto-hadronic scenario has also been suggested 
for the PWN in the SNR G54. 1+0.3 by Li et al. (2010). This 
source exhibits similar characteristics as GO.9+0.1 considering 
the characteristic age and spin-down luminosity of the pulsar. 
Finally it has to be noted that the IR and starlight seed-photon 
fields are only an estimate that does not contain local variations 
and could furthermore differ from the actual ones due to the un- 
certain distance to the PWN. 

7. Conclusion 

We have presented an extensive analysis of the nonthermal emis- 
sion of the PWN in the composite SNR GO.9+0. 1 . Unlike earlier 
publications, we were able to better resolve and determine spec- 
tral structures in the PWN thanks to an enhanced X-ray data set. 
Moreover, we performed the first spatially resolved modeling of 
this source in the X-ray regime and calculated the IC emission 
of the modeled part of the PWN. 

The hardness ratios of two annulus halves, separated by the 
presumed symmetry axis of the PWN, led us to conclude that 
we do not see any significant east-west deviation of the spectral 
properties of the source around the region of brightest emission. 
Thus we assume that PSR J1747-2809 is located in that region. 
To model the emission, we extracted spectra of four annulus- 
shaped regions centered on the putative pulsar position. We fixed 



the absorption column density to the value obtained from a fit of 
the whole modeled region as this approach significantly reduces 
the statistical error. The spectral index increases with larger dis- 
tance to the pulsar, whereas the surface brightness drops. 

These characteristics can be explained by assuming a lep- 
tonic outflow that suffers synchrotron and adiabatic energy 
losses as implemented in our model. We carried out the pa- 
rameter optimization using two different lepton injection spec- 
tra, namely a power law {Scenario I) and the injection spec- 
trum from Spitkovsky (2008) {Scenario II). Since the adopted 
model is only valid for the inner part of the PWN where the 
MHD limit is assumed to hold, we ignored the results of the 
fourth annulus. The evolution of the surface brightness can be 
reproduced for both injection spectra. While for Scenario I the 
slope of the modeled spectral index is far too low compared to 
the measured one, with Scenario II this characteristic can also 
be reproduced. Therefore we conclude that this injection spec- 
trum is better suited to explaining the nonthermal emission of 
this source, however, the energetics of the injection spectrum 
with the optimized parameters are shifted by about three orders 
of magnitude to higher values than in other works (see e.g. Slane 
201 1 ; Fang & Zhang 2010). The reason for this difference is the 
positive correlation of the spectral index of the injected lepton 
population and the one of the resulting synchrotron radiation. 
The power-law tail of the lepton spectrum (see Eq. 3) is sufficient 
for explaining the averaged spectral index of the whole PWN, as 
can be seen in these references. In contrast to that, our model 
is designed to explain the emission of several zones within the 
PWN, including the innermost one near the termination shock. 
Since the spectral index of that region is smaller, it is necessary 
to shift the spectrum to higher energies to get more of a contri- 
bution by the leptons of the relativistic Maxwellian. 

With the optimized model parameters it is possible to deduce 
some of the physical properties of the system. An example of this 
is the radial dependence of the magnetic field inside the PWN. 
Interestingly it increases with greater distance according to the 
model. However, an exact determination of the field strength is 
not possible since the model parameters are correlated, leading 
to multiple solutions of parameter sets with a low X 2 value. 

We calculated the IC emission of the modeled part of the 
PWN using the parameters optimized to reproduce the X-ray 
emission. The flux and spectral shape of the modeled radiation 
do not agree with the measured data. Several possibly concur- 
rent effects may offer an explanation. The unknown extent of 
GO.9+0.1 in VHE y-rays impedes an appropriate comparison of 
the measured and modeled data. That, in turn, makes it difficult 
to find out the relative contributions of the other effects, such 
as TeV emission from hadrons inside the PWN and older lep- 
tons that are not accounted for in our modeling. The problem of 
the unknown extent of GO.9+0.1 in VHE y-rays may be solved 
with future Imaging Atmospheric Cherenkov Telescopes like the 
Cherenkov Telescope Array (CTA). 
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